function y=Theorem3(u,mu,sigma,a,alpha,beta)
y=((2*mu/sigma^2-1)*(a-beta)^(2*mu/sigma^2-2)/u^(2*mu/sigma^2-1))*(((2*mu/sigma^2-1)*alpha^(-2*mu/sigma^2)/(a-beta)^(1-2*mu/sigma^2))/(1-(alpha/(a-beta))^(1-2*mu/sigma^2))+((1-2*mu/sigma^2)*(1-(alpha/a)^(1-2*mu/sigma^2))*alpha^(-2*mu/sigma^2)/((alpha/(a-beta))^(1-2*mu/sigma^2)))/(1-(alpha/(a-beta))^(1-2*mu/sigma^2))^2)+(1-(u/(a-beta))^(1-2*mu/sigma^2))*(-((2*mu/sigma^2-1)*alpha^(-2*mu/sigma^2)/(a-beta)^(1-2*mu/sigma^2))*((2*mu/sigma^2-1)*(a-beta)^(2*mu/sigma^2-2)/alpha^(2*mu/sigma^2-1))/(1-(alpha/(a-beta))^(1-2*mu/sigma^2))^2+(((1-(alpha/a)^(1-2*mu/sigma^2))*(1-2*mu/sigma^2)^2*alpha^(-2*mu/sigma^2)*(a-beta)^(2*mu/sigma^2-2))*(1-(alpha/(a-beta))^(1-2*mu/sigma^2))-2*((2*mu/sigma^2-1)*(a-beta)^(2*mu/sigma^2-2)/alpha^(2*mu/sigma^2-1))*((1-2*mu/sigma^2)*(1-(alpha/a)^(1-2*mu/sigma^2))*alpha^(-2*mu/sigma^2)/((alpha/(a-beta))^(1-2*mu/sigma^2))))/(1-(alpha/(a-beta))^(1-2*mu/sigma^2))^3);
% A=(1-(u/(a-beta))^(1-2*mu/sigma^2));
% B=(1-(alpha/(a-beta))^(1-2*mu/sigma^2));
% C=((1-2*mu/sigma^2)*(1-(alpha/a)^(1-2*mu/sigma^2))*alpha^(-2*mu/sigma^2)/((alpha/(a-beta))^(1-2*mu/sigma^2)));
% A'=((2*mu/sigma^2-1)*(a-beta)^(2*mu/sigma^2-2)/u^(2*mu/sigma^2-1));
% B'=((2*mu/sigma^2-1)*(a-beta)^(2*mu/sigma^2-2)/alpha^(2*mu/sigma^2-1));
% C'=((1-(alpha/a)^(1-2*mu/sigma^2))*(1-2*mu/sigma^2)^2*alpha^(-2*mu/sigma^2)*(a-beta)^(2*mu/sigma^2-2));
% k=((2*mu/sigma^2-1)*alpha^(-2*mu/sigma^2)/(a-beta)^(1-2*mu/sigma^2));
%y=A'*(k/B+C/B^2)+A*(-k*B'/B^2+(C'*B-2*B'*C)/B^3);